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ABSTRACT 



Aims. We present a unique multi-epoch infrared interferometric study of the oxygen-rich Mira variable RR Aql in comparison to 
radiative transfer models of the dust shell. We investigate flux and visibility spectra at 8 - 13yum with the aim of better understanding 
-^ , the pulsation mechanism and its connection to the dust condensation sequence and mass-loss process. 

'"J ■ Methods. We obtained 13 epochs of mid-infrared interferometry with the MIDI instrument at the VLTI between April 2004 and July 

2007, covering minimum to pre-maximum pulsation phases (0.45-0.85) within 4 cycles. The data are modeled with a radiative trans- 

(-H I fer model of the dust shell where the central stellar intensity profile is described by a series of dust-free dynamic model atmospheres 

O „ based on self-excited pulsation models. We examined two dust species, silicate and AI2O3 grains. We performed model simulations 

I . using variations in model phase and dust shell parameters to investigate the expected variability of our mid-infrared photometric and 

O ' interferometric data. 

Results. The observed visibility spectra do not show any indication of variations as a function of pulsation phase and cycle. The 

C/2 , observed photometry spectra may indicate intracycle and cycle-to-cycle variations at the level of 1-2 standard deviations. The pho- 

Cu ■ tometric and visibility spectra of RR Aql can be described well by the radiative transfer model of the dust shell thatuses a dynamic 

model atmosphere describing the central source. The best-fitting model for our average pulsation phase of Ov = 0.64 ±0.15 includes 

the dynamic model atmosphere M21n (rmodei = 2550 K) with a photospheric angular diameter of Sphoi = 7.6 + 0.6 mas, and a silicate 

►^ . dust shell with an optical depth of Ty = 2.8 ± 0.8, an inner radius of R,„ = 4.1 ± 0.7i?phot, and a power-law index of the density 

r^v ' distribution of p = 2.6 ± 0.3. The addition of an AI2O3 dust shell did not improve the model fit. However, our model simulations 

^— s ' indicate that the presence of an inner AI2O3 dust shell with lower optical depth than for the silicate dust shell can not be excluded. 

r^ , The photospheric angular diameter corresponds to a radius of iJphot = 520*'^^Ro and an effective temperature of T^if ~ 2420 + 200 K. 

^— V , Our modeling simulations confirm that significant intracycle and cycle-to-cycle visibility variations are not expected for RR Aql at 

• ' mid-infrared wavelengths within our uncertainties. 

P^ ' Conclusions. We conclude that our RR Aql data can be described by a pulsating atmosphere surrounded by a silicate dust shell. The 

(«_^ I effects of the pulsation on the mid-infrared flux and visibility values are expected to be less than about 25% and 20%, respectively, 

and are too low to be detected within our measurement uncertainties. 
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;h 1 ■ Introduction very cool and dense outer part of the pulsating atmosphere. Dust 
^ ] created by carbon-rich stars owing to its high opacity absorbs 
Asymptotic giant branch (AGB) stars are low-to-intermediate radiation from the star, and the wind is believed to be driven by 
mass stars at the end of their stellar evolution. These stars, in- radiation pressure on the dust particles that drag gas along via 
eluding Mira type stars, exhibit many complex processes such as collisions ("Hofner 20081 
shock fronts propagating through the stellar atmosphere, large- 
amplitude pulsation, and molecule and dust formation leading to Despite remarkable progress in theoretical and observational 
strong mass loss via a dense and dusty outflow from an extended studies, it still remains unclear whether the dust in oxygen-rich 
stellar atmosphere (Andersen et al. 2003). The mass loss rates envelo pes is opaque enough to d rive the observed massive out- 
can reach up to 10 '^Mgyr ' (Matsuura et al. 2009) with expan- flows ([Woitkell2006t lHofneill2008 1. Because of the wind, all the 
sion velocities of 5-30 km s ' (Hofner & Anderseni2007.) . Dust matter around the core of the star is eventually returned to the 
grains are created by condensation from the gas phase in the interstellar medium, and the star evolves toward the planetary 
nebula phase, leaving the core as a white dwarf. Generally many 



* Based on observations made with Very Large Telescope aspects of the physics of AGB stars remain poorly understood. 
Interferometer (VLTI) at the Paranal Observatory under program IDs such as the detailed atmospheric stratification and composition 
073.D-0711, 075.D-0097, 077.D-0630, and 079.D-0172. of the stars, or the role of stellar pulsation and its connection 
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to the dust formation and massive outflow. Improving our un- 
derstanding of the physical processes leading to the substantial 
niass loss is important, as AGB stars play a crucial role in the 
chemical enrichment and evolution of galaxies by returning gas 
and dust to the interstellar medium (ISM). 

Thanks to their large diameters and high luminosities, Mira 
variables are ideal targets for high angular resolution obser- 
vations. Near-infrared interferometric (NIR) techniques pro- 
vide detailed information regarding the conditions near the 
continuum-forming photosphere such as effective temperature, 
center-to-limb intensity variations, the stellar photospheric di- 
ame ter, and its dependence on wavelength and pulsation phase 
(e.g. lHaniffetal.ll 19951 iPerrin et al.lll999l; Ivan Belle et a l.lll996t 
Young et al. 2000a; Thompson et al. 2002a b; Ouirrenba ch et alJ 
19921 " iFedele et al. 2005; Ohnakai i2004t iMillan-Gabet et alj 
200I [Woodruff etal. 20081). New theoretical self-excited dy- 



Stellar cycle/phase 



namic model atmospheres of o xygen-rich stars have been cre- 
ated and successf ully appli ed ( Tei et al. 20031; iHofma nn et al] 
ll998Hlreland"etai . 2004b a. l2008HWittkowski e t al. 200l^ 

Observations at mid-infrared wavelengths are well-suited to 
studying the molecular shells and the dust formation zone of 
evolved stars. Mid-infrared interferometry is sensitive to the 
chemical composition and geometry of dust shells, their tem- 
pera ture, inner radii, radial distribution, and the mass loss rate 
(e.g.'Danchi et al.'l994';'Monnier et al.ll997l l2000';'Wein er et alJ 
12006) . Lopez et al. ( 1997) presented long-term observations at 
1 1 fivn of o Ceti obtained with the Infrared Spatial Interferometer 
(ISI). The observed visibilities change from one epoch to the 
next and are not consistent with simple heating or cooling of the 
dust with change in luminosity as a function of stellar phase, 
but rather with large temporal variations in the density of the 
dust shell. The data were compared to axially symmetric ra- 
diative transfer models and suggest inhomogeneities or clumps. 
iTevousia net al.l ([2004) has studied the spatial distribution of 
dust around four late-type stars with ISI at 11.15yum, and find 
that the visibility curves change with the pulsation phase of 
the star. The dust grains were modeled as a mixture of sili- 
cates and graphite. The results suggest that the dust shells ap- 
pear to be closer to the star at minimum pulsation phase, and 
farther away at maximum phase, which is also demonstrated by 
IWittkowski et al.1 (l2007l) . lOhnaka et alj (l2007h reported on tem- 
poral visibility variations of the carbon-rich Mira variable V Oph 
with the instruments VINCI and MIDI at the VLTI. This tempo- 
ral variation of the N-band angular sizes is largely governed by 
the variations in the opacity and the geometrical extension of the 
molecular layers (C2H2 and HCN) and the dust shell (amC + 
SiC). 

Additional information to the near- and mid-infrared obser- 
vations can be obtained with Very Long Baseline Array (VLBA) 
observations, which allow one to study the properties of the cir- 
cumstellar environment of evolved stars using the maser radia- 
tion emitted by some molecules, most commonly SiO, H2O, and 
OH. By observing the maser radiation at tens to several hundred 
AU from the star, it is possible t o determine a ccurate proper mo- 
tions and parallaxes ( Vlemming s et al. 2003h . total intensity and 
linear polarizations (Cotton et al. '2008"), and the structure of the 
environment of the stars (Diamond et al. 1994). 

This work is part of the ongoing project of concurrent mul- 
tiwavelength observations using optical long baseline interfer- 
ometry (AMBER and MIDI at the VLTI) combined with ra- 
dio interferometry (VLBA). We aim to investigate layers at 
different depths of the atmosphere and circumstellar envelope 
(CSE) of the stars. Previous results from this project of joint 
VLTI/VLBA observations for the Mira star S Ori can be found 
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Fig. 1. Visual light curve of RR Aql based on data from the 
AAVSO and AFOEV databases as a function of Julian Date and 
stellar cycle/phase. The green arrows indicate the dates of our 
VLTI/MIDI observations. The blue arrows indicate the dates of 
our VLBA observations. Here, the full arrows denote observa- 
tions of SiO maser emission and simple arrows observations of 
H2O maser emission. 



in iBoboltz & Wittkowskil (l2005h and IWittkowski eTal] (l2007h . 
Here we present results of a long-term VLTI/MIDI monitor- 
ing of the oxygen-rich Mira variable RR Aql. In addition to 
the VLTI/MIDI observation presented in this paper, we obtained 
coordinated multi-epoch VLBA observations of SiO and H2O 
masers towards RR Aql, which will be presented in a subsequent 
paper. 

2. Characteristics of RR Aql 

RR Aql is an oxygen- r ich M ira variable with spectral type 
M6e-M9 (ISamus et al.1 |2004|) . RR Aql shows a sti'ong 
silicate emission feature in i ts mid-infrared spectrum 
(JLorenz-Martins & Pompeial l2000l) . In addition, it has rela- 
tively s trong SiO, H2O, and OH maser emission dBenson et al.l 
Il990h . IVlemmings & van Langeveldd ( l2007l) observed OH 
masers toward RR Aql at five epochs with the VLBA. Based 
on these observations, the distance to RR Aql is estimated to 
D = 633!^i4pc (HIPPARCOS distance D = 540 pc). Van Belle 
et al. (2002) obtained a /T-band (A = 2.2 yum, AA = 0.4 yum) 
angular size with the Infrared Optical Telescope Array (IOTA). 
The uniform disk diameter (UD) of 0ud = 10.73 + 0.66 mas at 
phase O = 0.48, an effective temperature T^ff = 2127 ± 111 K, 
and a bolometric fl ux /boi = 78.4 ±11.8 10"^ ergs cm^^s"'. 
'Mivata et al.l (l2000 i) studied the dust around the star spectro- 
scopicall y and obt ained Fdust/^star at lOjum of 1.49 + 0.02. 
Ragland et al.l (12006 ) measured the closure phase with IOTA in 
the //-band, and classified RR Aql as a target with no detectable 
asymmetries. The IRAS flux at 12jum is 332 Jy. The light curve 
in the V band varies from ~8 mag at maximum light to ~ 14 mag 
at minimum light. IWhitelock et al.l (2000) shows that the K 
magnitude of RR Aql varies from ~0.0 mag at maximum light to 
~ 1 .0 mag at mini mum light. RR Aql i s pulsating with a period of 
P = 394.78 days dSamus et al.ll2004 and the J ulian Date of flie 
last maximum brightness is Tq - 2452875.4 ((Poimans ki et al.l 
2005). The period of pulsation visually corresponds to recent 
values from the AAVSCQ and AFOEVQ databases. Figure [T| 



' http://www.aavso.org 

" http://cdsweb.u-strasbg.fr/afoev 
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Table 1. VLTI/MIDI observation of RR Aql. 



Ep. 


date 
DD/MMrfYVY 


Time 
[UTC] 


JD 


*,.,. 


Conf. 


B 


Disp. 
Elem. 


BC 


Br 
[m] 


PA. 
[deg] 


Seeing 
["] 


To 

[msec] 


A 


09/04/2004 


07:26 


2453105 


0.58 


U2-U3 


Aim 


Prism 


HS 


33.57 


3.21 


0.45 


7.4 


A 


10/04/2004 


09:05 


2453106 


0.58 


U2-U3 


Aim 


Prism 


HS 


37.17 


26.09 


0.76 


7.8 


A 




09:52 


2453106 


0.58 


U2-U3 


Aim 


Prism 


HS 


39.86 


33.48 


0.64 


8.2 


B 


09/07/2004 


08:25 


2453196 


0.81 


U2-U3 


Aim 


Prism 


HS 


45.23 


44.83 


0.47 


4.8 


C 


28/07/2004 


08:09 


2453215 


0.86 


U2-U3 


Aim 


Prism 


HS 


42.54 


41.57 


1.36 


1.0 


C 


29/07/2004 


02:05 


2453216 


0.86 


U2-U3 


Aim 


Prism 


HS 


37.83 


28.16 


1.66 


1.0 


C 




03:33 


2453216 


0.86 


U2-U3 


Aim 


Prism 


HS 


42.84 


39.49 


0.75 


2.3 


C 




04:55 


2453216 


0.86 


U2-U3 


Aim 


Prism 


HS 


45.99 


44.67 


1.04 


1.6 


C 




06:33 


2453216 


0.86 


U2-U3 


Aim 


Prism 


HS 


46.16 


45.71 


3.38 


0.5 


C 




07:04 


2453216 


0.86 


U2-U3 


Aim 


Prism 


HS 


45.32 


44.92 


2.68 


0.7 


C 


01/08/2004 


00:12 


2453219 


0.87 


U2-U3 


Aim 


Prism 


HS 


33.76 


6.86 


0.73 


3.1 


D 


18/04/2005 


10:30 


2453479 


1.53 


U2-U4 


89m 


Prism 


HS 


88.37 


81.72 


0.76 


3.5 


D 


19/04/2005 


07:53 


2453480 


1.53 


U2-U4 


89m 


Prism 


HS 


61.36 


76.41 


0.83 


2.3 


E 


20/07/2005 


03:51 


2453572 


1.76 


U1-U4 


130m 


Prism 


HS 


119.63 


59.78 


0.71 


1.9 


E 




06:16 


2453572 


1.76 


U1-U4 


130m 


Prism 


HS 


128.64 


63.28 


0.81 


1.6 


E 


22/07/2005 


05:10 


2453574 


1.77 


U2-U3 


47m 


Prism 


HS 


45.65 


44.12 


1.00 


1.4 


E 




07:56 


2453574 


1.77 


U2-U3 


47m 


Prism 


HS 


44.44 


43.96 


1.11 


1.2 


F 


18/04/2006 


09:38 


2453844 


2.45 


DO-GO 


32m 


Prism 


HS 


28.99 


70.01 


1.08 


2.1 


F 


19/04/2006 


08:32 


2453845 


2.46 


DO-GO 


32m 


Prism 


HS 


24.74 


65.81 


0.80 


6.0 


F 




09:21 


2453845 


2.46 


DO-GO 


32m 


Prism 


HS 


28.21 


69.32 


0.78 


6.3 


G 


21/05/2006 


05:53 


2453877 


2.54 


EO-GO 


16m 


Prism 


HS 


10.94 


62.10 


0.56 


4.0 


G 


25/05/2006 


06:35 


2453881 


2.55 


AO-GO 


64m 


Prism 


HS 


53.17 


67.77 


0.55 


5.1 


G 




08:46 


2453881 


2.55 


AO-GO 


64m 


Prism 


HS 


63.93 


72.67 


0.44 


6.7 


H 


18/06/2006 


03:58 


2453905 


2.61 


DO-GO 


32m 


Prism 


HS 


21.37 


61.34 


1.29 


1.5 


H 




04:41 


2453905 


2.61 


DO-GO 


32m 


Prism 


HS 


25.09 


66.21 


1.27 


1.6 


H 


18/06/2006 


05:39 


2453905 


2.61 


DO-GO 


32m 


Prism 


HS 


29.05 


70.06 


1.00 


2.1 


H 


20/06/2006 


05:49 


2453907 


2.61 


AO-GO 


64m 


Prism 


HS 


59.89 


70.81 


1.05 


3.9 


H 




06:33 


2453907 


2.61 


AO-GO 


64m 


Prism 


HS 


63.00 


72.15 


0.88 


4.5 


H 


21/06/2006 


07:08 


2453908 


2.62 


EO-GO 


16m 


Prism 


HS 


16.00 


72.76 


1.01 


2.2 


H 




09:36 


2453908 


2.62 


EO-GO 


16m 


Prism 


HS 


13.21 


70.37 


0.77 


2.9 


H 


23/06/2006 


03:20 


2453910 


2.62 


EO-GO 


16m 


Prism 


HS 


9.84 


58.42 


0.76 


3.9 


H 




04:07 


2453910 


2.62 


EO-GO 


16m 


Prism 


HS 


11.97 


64.86 


0.67 


4.4 




08/08/2006 


05:12 


2453956 


2.74 


AO-GO 


64m 


Prism 


HS 


61.41 


72.62 


1.02 


2.3 






05:57 


2453956 


2.74 


AO-GO 


64m 


Prism 


HS 


56.88 


71.58 


1.10 


2.2 




09/08/2006 


06:39 


2453957 


2.74 


AO-GO 


64m 


Prism 


HS 


50.33 


69.49 


1.47 


1.5 




10/08/2006 


03:07 


2453958 


2.74 


DO-GO 


32m 


Prism 


HS 


31.37 


72.03 


1.43 


2.6 






04:31 


2453958 


2.74 


DO-GO 


32m 


Prism 


HS 


31.65 


72.90 


1.67 


2.5 






05:45 


2453958 


2.74 


DO-GO 


32m 


Prism 


HS 


28.72 


71.74 


1.43 


3.2 




11/08/2006 


02:44 


2453959 


2.74 


DO-GO 


32m 


Prism 


HS 


30.76 


71.51 


1.43 


2.3 




13/08/2006 


04:05 


2453961 


2.75 


EO-GO 


16m 


Prism 


HS 


15.94 


72.90 


0.91 


1.9 






04:56 


2453961 


2.75 


EO-GO 


16m 


Prism 


HS 


15.26 


72.56 


0.63 


2.7 






05:52 


2453961 


2.75 


EO-GO 


16m 


Prism 


HS 


13.76 


71.07 


0.60 


2.8 




14/09/2006 


01:01 


2453993 


2.83 


EO-GO 


16m 


Prism 


HS 


15.82 


72.29 


0.85 


2.3 




16/09/2006 


02:15 


2453995 


2.84 


AO-GO 


64m 


Prism 


HS 


62.86 


72.85 


0.95 


1.6 






03:13 


2453995 


2.84 


AO-GO 


64m 


Prism 


HS 


58.12 


71.90 


1.68 


0.9 


K 


09/05/2007 


08:53 


2454230 


3.43 


DO-HO 


64m 


Prism 


HS 


61.43 


71.46 


1.12 


1.3 


L 


22/06/2007 


03:50 


2454274 


3.54 


GO-HO 


32m 


Prism 


HS 


21.96 


62.20 


1.17 


1.0 


L 




05:29 


2454274 


3.54 


GO-HO 


32m 


Prism 


HS 


29.30 


70.26 


1.21 


0.9 


M 


03/07/2007 


07:25 


2454285 


3.57 


EO-GO 


16m 


Prism 


HS 


15.51 


72.73 


1.07 


2.8 


M 




07:37 


2454285 


3.57 


EO-GO 


16m 


Prism 


HS 


15.31 


72.59 


0.75 


3.9 


M 


04/07/2007 


04:42 


2454286 


3.57 


GO-HO 


32m 


Prism 


HS 


29.28 


70.24 


1.59 


1.0 


M 




06:30 


2454286 


3.57 


GO-HO 


32m 


Prism 


HS 


31.98 


72.83 


1.51 


1.1 



Notes. The table lists the epoch, the date, the time, the Julian Date (JD), the visual pulsation phase O,,,,, the baseline configuration, the ground 
length of the configuration, the dispersive element, the beam combiner BC, the projected baseline length Bp, the position angle on the sky P.A. 
(deg. east of north), the DIMM seeing (at 500 nm), and the coherence time Tq (at 500 nm). 



shows the visual light curve of RR Aql based on values from the 
AAVSO and AFOEV databases as a function of Julian day and 
stellar phase. 



3. VLTI/MIDI Observations and data reduction 

We obtained 57 spectrally-dispersed mid-infrared interferomet- 
ric observations of RR Aql with the VLTI/MIDI instrument be- 
tween Apr 9, 2004 and Jul 28, 2007. MIDI - th e mid-infrared 
interferometric instrument dLeinert et al.l |2004|) combines the 
beams from two telescopes of the VLTI dGlindemann et al.) 
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Fig. 2. Calibrated RR Aql MIDI flux spectra as a function of 
wavelength. For clarity, the error bars are omitted in the plot. 



l2003h and provides spectrally resolved visibilities in the lOyum 
window (A^-band, 8-13/im). To obtain dispersed photometric 
and interferometric signals, we used the PRISM as a dispersive 
element with a spectral resolution R = AA/A ~ 30. The beams 
were combined in HighJSens mode (HS). In this mode the pho- 
tometric signal is observed after the interferometric signal. 

The details of the observations and the instrumental settings 
are summarized in Table [T] All observations were executed in 
service mode using either the unit telescopes (UTs, 8.2m) or the 
auxiliary telescopes (ATs, 1.8m). 

We merged the MIDI data into 15 epochs, with a maximum 
time-lag between individual observations of five days for each 
epoch (1.3% of the pulsation period). For technical problems we 
had to eliminate two epochs. Fig.[T]shows the final 13 epochs in 
comparison to the light curve. The uncertainty in the allocation 
of the visual phase to our observations was estimated to ~0. 1 . 
Our long-term VLBA monitoring of SiO and H2O maser emis- 
sion toward RR Aql, will be presented in a forthcoming paper. 

We used the MIA+EWS software package, version 1.6 
(Jaffe, Koehler, et aQ} for the MIDI data reduction. This package 
includes two different methods, an incoherent method (MPIA 
software package MIA) that analyzes the powerspectrum of the 
observed fringe signal and a coherent integration method (EWS), 
which first compensates for optical path differences, including 
both instrumental and atmospheric delays in each scan, and then 
coherently adds the fringes. We applied both methods to inde- 
pendently verify the data reduction results. The detector masks 
were calculated by the procedure of MIA, and were used for 
both the MIA and the EWS analysis. The obtained data reduc- 
tion values correspond to each other, and we chose to use in the 
following the results derived from the EWS analysis, which offer 
error estimations. 

To account for instrumental visibility losses and to deter- 
mine the absolute flux values, calibrator stars with known flux 
and diameter values were observed immediately before or af- 
ter the science target. Our main calibrators were HD 169916 
(period P73, P75), HD 146051 (P77), and HD 177716 (P79). 
Calibrated science target visibility spectra were calculated us- 
ing the instrumental transfer function derived from all calibrator 
data sets taken during the same night with the same baseline and 
instrumental mode as our scientific target. The number of avail- 
able transfer function measurements depends on the number of 
calibrator stars observed per specific night including those cal- 
ibrators observed by other programs. The errors of the transfer 
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Fig. 3. RR Aql flux spectrum as a function of wavelength from 
2.4//m to 40 //m (top). The lines represent the flux spectra from 
ISO-SWS observations covering wavelengths from 2.5 yum to 
~240/im (solid red thick line), IRAS-LRS observations cover- 
ing ~ 7.7 jum to ~ 23jum (dashed thin black line), and the mean 
of our A^-band MIDI measurements (solid thin green line). Here, 
the vertical bars span the maximum and minimum values mea- 
sured. The diamonds denote 2MASS measurements at 1.25//m, 
1.65/im, and 2.2 fim. The solid black line indicates our atmo- 
sphere and dust shell model as explained in Sect. 15. 21 The bottom 
plot shows an enlarged segment of the plot in the MIDI wave- 
length range of 8-13jum. 

functions are given by the standard deviation of all transfer func- 
tion measurements per night. To estimate the uncertainty of the 
transfer function for nights when only one calibrator was avail- 
able, we used typical values based on nights when many calibra- 
tor stars were observed. The final errors on the observed visibili- 
ties are mostly systematic, and include the error of the coherence 
factor of the science target and the calibrators, the adopted di- 
ameter errors, and the standard deviation of the transfer function 
over the night. 

The photometric spectrum was calibrated with one or two 
calibrators, which were observed close on the sky and in time 
compared to our science target. For most calibrator stars, abso- 
lutely calibrated spectra are available in Cohen et al. ( 1999). For 
those calibrators where the absolutely calibrated spectrum was 
not directly available, we instead used a spectrum of a calibra- 
tor with a similar spectral type and similar effective temperature 
(see the instrument consortium's catalog^). The spectra of such 
calibrators were scaled with the IRAS flux at 12yum to the level 
of our calibrator. In addition, we verified that our synthetic spec- 
tra obtained by this procedure are valid by scaling known spectra 



http://www.strw.leidenuniv.nl/~nevec/MIDI 



http://www.ster.kuleuven.ac.be/~tijl/MIDI_calibration/mcc.txt 
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Fig. 4. Calibrated RR Aql MIDI visibility amplitudes as a function of wavelength, (top left) Observations executed at projected 
baselines (Bp) < 20 m, (top right) 20 m < Bj, < 35 m, (bottom left) 35 m < Bp < 50 m, and (bottom right) Bp > 50 m. 



of two Cohen calibrators. In a few cases, when the atmospheric 
absorption was strongly affecting the spectra around 9.5/im, we 
used another similar calibrator instead of the main calibrator ob- 
served in the same night with the same level of flux. The ambi- 
ent conditions for all the observations were carefully checked. If 
we were facing problems with clouds, constraints due to wind, 
significant differences in seeing, humidity, coherence time, and 
airmass between the science target and the corresponding cali- 
brator, the photometry was omitted from the analysis. 



4. MIDI results 

We observed RR Aql at different phases over four pulsation cy- 
cles in order to monitor the mid-infrared photometry and visibil- 
ity spectra. Here, we present and discuss the general properties 
of the data, followed by an analysis of their variability as a func- 
tion of phase and cycle in Sects. 143] and l4. 21 

Figure|2]shows all obtained calibrated photometry spectra as 
a function of wavelength. The MIDI flux measurements show a 
consistent shape exhibiting an increase in the flux from ~ 100- 
200 Jy at 8jum to a maximum near 9.8 jum of ~200-400Jy, and a 
decrease towards I3fj.m, where the flux values again reach val- 
ues of ~ 1 00-200 Jy. The level of the flux spectrum differs for in- 
dividual measurements with a spread of ~ 100-200 Jy. Figure[3] 
shows a comparison of the mean of our MIDI flux measure- 
ments to measurements obtained with the ISO and IRAS in- 
struments. The shape of the flux curve is consistent among the 
MIDI, ISO, and IRAS measurements. The level of the IRAS flux 
is within the range of our MIDI measurements. The level of the 
ISO flux is higher, which can most likely be explained by the 
post-maximum phase of 0.16 of the ISO observations (1997-05- 



03) compared to our minimum to pre-maximum phases. The dif- 
ference between the MIDI and ISO flux level may also indicate 
a flux variation over different cycles. Figure [3] also includes a 
model description of our MIDI data, which is explained below 
in Sect. EH 

Figure |4] shows all obtained calibrated visibility spectra, 
which are combined into four groups of different projected base- 
Hne lengths (Bp) of Bp < 20 m, 20 m < B^ < 35 m, 35 m < Bp 
< 50 m, and Bp > 50 m. The visibility curves show a signifi- 
cant wavelength dependence with a steep decrease from 8 yum to 
~ 9.5 jum and a slow increase in the 9.5 //m to 13jum range. The 
shape and absolute scale of the visibility function depends on the 
length of the projected baseline where a longer baseline results 
in a lower and flatter curve. 

As a first basic interpretation of the interferometric data, we 
computed the corresponding uniform disk (UD) diameter and the 
Gaussian FWHM for each data set and spectral channel. This 
leads to a rough estimate of the characteristic size of the tar- 
get at each wavelength. However, it should be mentioned that 
the true intensity distribution across the stellar disk is expected 
to be more complex than can be described by these elemen- 
tary models. Figure |5] shows for the example of epoch A, and 
epoch I the flux, the visibility amplitude, the corresponding UD 
diameter, and the corresponding Gaussian FWHM diameter as 
a function of wavelength. The shape of the UD diameter and 
Gaussian FWHM functions show a steep increase from 8 yum to 
~9.5/im by a factor of 2. A plateau appears between ~9.5/im 
and 11.5 fj.m, and the values are nearly constant from ~ 1 1.5jum 
to 13//m. 

The shape of the 8-13/im flux spectrum with a max- 
imum near 9.8 //m is known to be a characteristic sil- 
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Fig. 6. Calibrated MIDI visibility amplitudes for different pulsa- 
tion phases within the same cycle to investigate intracycle visi- 
bility variations. Each line represents a different pulsation phase 
within the same cycle and is computed as an average of data ob- 
tained at the respective phase (+ 0.15) and observed at similar 
projected baseline length (Bp + 10%) and position angle (P. A. 
+ 10%). The top panel shows the example of pulsation phases 
0.45, 0.61, and 0.74 of cycle 2 observed with a projected base- 
line length of ~ 29 m and a position angle of ~70deg. The bottom 
panel shows the example of pulsation phases 0.62, 0.75, and 0.83 
of cycle 2 observed with a projected baseline length of ~ 16 m 
and a position angle of ~72deg. The error bars are computed as 
the standard deviation of the averaged visibilities. 



icate emission feature (e.g. iLittle-M arenin & Littla Il990t 
iLorenz-Martins & Pompeia 12000) . The corresponding drop in 
the visibiUty function between 8jum and 9.5 jum and increase in 
the UD diameter and Gaussian FWHM has also been shown to 
be a typical signature of a circu mstellar dust shel l that is domi- 
nate d by silicate dust (cf., e.g. Drieb e et al]|2008l; lOhnaka et al.l 
12008 ). The opposite trend toward a broad spectrum and visibil- 
ity function in the wavelength range 8-13 ^m was interpreted as 
a dust shell consisting of AI2O3 dust (Witt kowski et al.ll2007h or 
silicate and AI2O3 dust (.Qhnaka et al.,2005lK 

4. 1 . Visibility monitoring 

We obtained a rich sample of MIDI data on RR Aql, which 
covers a total of 4 pulsation cycles and pulsation phases be- 
tween 0.45 and 0.85, i.e. minimum to pre-maximum phases (see 
Table [1] and Fig. [TJ. For many different pulsation cycles and 
phases, we obtained data at similar projected baseline lengths 
and position angles. This gives us a unique opportunity to mean- 
ingfully compare interferometric data obtained at different pul- 



Fig. 7. As for Fig. |6] but for the same pulsation phase in con- 
secutive pulsation cycles to investigate cycle-to-cycle visibility 
variations. The top panel shows the example of pulsation phase 
~ 0.5 in three consecutive cycles observed with a projected base- 
line length of ~ 61m and a position angle ~75deg. The bottom 
panel shows the example of pulsation phase 0.5 in 2 consecutive 
cycles observed with a projected baseline length of ~ 29 m and 
a position angle ~70deg. 



sation phases and cycles. Since the visibility depends on the 
probed point on the uv plane, and thus on the projected base- 
line length and position angle, visibility data at different pulsa- 
tion phases can be directly compared only if they were obtained 
at the same, or very similar, point on the uv plane. For this rea- 
son, we combined individual observations into groups of similar 
pulsation phases (<l>vi.s ± 0.15), projected baseline lengths (Bp + 
10 %), and position angles (PA. + 10 %). The data within each 
group were averaged. The uncertainty of the averaged visibility 
curves was estimated as the standard deviation of the averaged 
visibilities. 



4.1.1. Intracycle visibility monitoring 

To investigate intracycle visibility variations we compared data 
observed at different pulsation phases within the same cycle. 
Figure |6] shows two examples of calibrated visibility curves, 
where each line in the plot represents an average of a group of 
visibility data, as described above, (i.e. with similar Bp, PA., and 
Ovis). The top panel shows an example of observations within the 
second pulsation cycle at phases 2.45, 2.61, and 2.74 obtained 
with Bp ~ 30 m and PA. ~ 70°. The bottom panel shows an ex- 



ample of observations at phases 2.62, 2.75, and 2.83 obtained 
with Bp ~ 16 m and P. A. ~ 72°. The data do not show any evi- 
dence of intracycle visibility variations within the probed range 
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Fig. 8. Calibrated MIDI flux spectrum for different pulsation 
phases within the same cycle to investigate intracycle photom- 
etry variations. Each line represents a different pulsation phase 
within the same cycle and is computed as an average of data ob- 
tained at the respective phase (+ 0.15). The top panel shows the 
example of pulsation phases 0.58 and 0.85 of cycle 0. The bot- 
tom panel shows the examples of phases 0.53 and 0.76 of cycle 
1 . The error bars are computed as the standard deviation of the 
averaged photometry curves. The gray shades denote zones that 
are affected by atmospheric absorption. 



of pulsation phases (~ 0.45 
curacies of about 5-20%. 



0.85) and within our visibiUty ac- 



4.1.2. Cycle-to-cycle visibility monitoring 

We compared data observed at similar pulsation phases of dif- 
ferent consecutive cycles to investigate cycle-to-cycle visibility 
variations. Figure Q shows two of these examples computed in 
the same way as shown in Fig.|6]of Sect. 14.1. Tl but (top panel) 
for minimum phases of 1.53, 2.59, and 3.43 in three consecutive 
cycles obtained with a projected baseline length of ~ 62 m and 
P.A. ~ 73°, and (bottom panel) for minimum phases 2.43, and 
3.55 of two consecutive cycles obtained with a projected base- 
line length of ~ 29 m, and P.A. ~ 70°. Here, the phases of the 
different cycles differ by up to ~ 20%. However, we showed in 
Sect. 14. l.T] that there is no evidence of intracycle visibility varia- 
tions, so that the chosen phases can be compared well for differ- 
ent cycles. As a result, Fig.|2]shows that our data do not exhibit 
any significant cycle-to-cycle visibility variation for minimum 
phases over two or three cycles within our visibility accuracies 
of ~ 5-20%. 
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Fig. 9. As for Fig. [8] but for the same pulsation phase in consec- 
utive pulsation cycles to investigate cycle-to-cycle photometry 
variations. The top panel shows the example of phase ~ 0.75 in 
cycles 1 and 2. The bottom panel shows the example of phase ~ 
0.85 of cycles and 2. 



4.1.3. Deviations from circular symmetry 

Most position angles varied around 40° or 70°. Only in two 
cases is the P.A. significantly different with values of 3° and 7°. 
However, in both cases the visibility data quality is poor. The 
differential phases are close to zero within ~ 10-20deg, which 
is close to the calibration unce rtainties that c an be reached for 
MIDI differential phases (e.g. lOhnakaet al.1 12008). Therefore 
these observations do not allow us to make any reliable conclu- 
sion about the presence or absence of an asymmetric intensity 
distribution. 

Summarizing, our RR Aql MIDI data do not show any evi- 
dence of significant intracycle or cycle-to-cycle visibility varia- 
tions within the examined phase coverage of ~ 0.45-0.85 in four 
consecutive cycles (see Table[T]for details of the phase coverage) 
and within our visibility accuracies of ~ 5-20%. The obtained 
mid-infrared interferometric observations imply either that the 
mid-infrared sizes of molecular and dusty layers of RR Aql 
do not significantly vary within our phase coverage or that the 
conducted observations are not sensitive enough to detect such 
variations. In addition, the good agreement of visibility data ob- 
tained at the same projected baseline lengths and position angles 
for several years also confirms the good data quality and credi- 
bility of the data reduction procedure. 
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4.2. N-band flux monitoring 

To investigate the intracycle and cycle-to-cycle variability of 
our 8-1 3 /vm flux spectra, we combined the data into groups of 
data obtained at similar pulsation phases (Ovis ± 0.15). The data 
within each group were averaged. The uncertainty of the aver- 
aged photometry curves was estimated as the standard deviation 
of the averaged values. 



4.2.1. Intracycle photometry monitoring 

Figure [8] shows two examples of a comparison of the MIDI cal- 
ibrated flux spectrum at different pulsation phases of the same 
cycle. The top panel shows a comparison between phases 0.58 
and 0.85 of cycle 0, and the bottom panel shows a compari- 
son between phases 1.53 and 1.76 of cycle 1. Both examples 
compare data taken at post-maximum and pre-minimum phases, 
i.e. with the largest separation in phase that is available in our 
dataset. The RR Aql flux values within the A^-band are larger at 
the pre-maximum phases than at the post-minimum phases, cor- 
responding to the light curve of the V-magnitude. The difference 
is most pronounced toward the silicate emission feature around 
9.8 yum and smaller towards the edges of the MIDI bandpass at 
8.0/im and 13.0/im. This can most likely be understood as a 
consequence of the silicate absorption coefficient being largest 
at short visual wavelengths, so that the silicate temperature is 
particularly sensitive to the large visual flux variations of the 
central star. The differences correspond to up to 20%-35%, or 
about one or two standard deviations. 



4.2.2. Cycle-to-cycle photometry monitoring 

Figure |9] shows a comparison of data observed at the same 
phase of consecutive pulsation cycles. The top panel shows 
a comparison between phases 1.76 and 2.74, i.e. phase 0.75 
of cycles 1 and 2. The bottom panel shows a comparison 
between phases 0.85 and 2.84, i.e phase 0.85 of cycles and 
2. In both examples, the A^-band fluxes are lower in cycle 
2 than for the same phase of (a) cycle I and (b) cycle by 
up to ~ lOOJy, corresponding to about 30% or one to two 
standard deviations. As for the intracycle photometry variations 
discussed in the previous paragraph, the differences are most 
pronounced toward the silicate emission feature at 9.8 jum. This 
might indicate an irregularity in the variability cycle in the 
A^-band. In the optical band the visual light curve clearly shows 
that the pulsation cycles are not perfectly symmetric (see Fig.[T]). 

Summarizing, our data exhibit a l-2cr signature of intracy- 
cle, as well as cycle-to-cycle, flux variations at wavelengths of 
8-1 3 //m, which are most pronounced toward the silicate emis- 
sion featu re at 9.8 Atm. This ind ication is consistent with obser- 
vations bv lMonnier et alj (Il998h . who report temporal variations 
in the mid-infrared spectra of late type stars. In particular, stars 
with a strong emission feature around 9.8 jum showed evident 
changes in the spectral profile. Considera ble phase variation of 
the A^-band spectra can also b e found in .Alvarez & Pled (Il998l) 
and lGuha Nivogi et alJ (l20Ilh . 

5. Modeling of our MIDI data 

Self-consistent models describing the dynamic atmosphere of 
the central source and the dust shell are still very rare. This is 
particularly true for oxygen-rich stars. Nevertheless, some ad- 
vance in this domain has already been successfully achieved 



(llreland & Scholzl2006HH"ofner & Andersenl2007Hlreland et al.1 



120081 : lHofneill2008 l). Here, we use an approach of an ad hoc 
radiative transfer modeling of the dust shell where the central 
stellar source is described by readily available and established 
dust-free dynami c model atmosphere series, as introduced by 
IWittkowski eTai] (|2007). 

The dust shell surrounding the central star is modele d by the 
Monte Carlo radiative transfer code mcsimjnpi ( Ohnak a et alj 
|2006). The radiative transfer model requires an assumption on 
the spectral energy distribution (SED) of the central stellar 
source. Since we expect the photosphere and molecular layers 
of RR Aql to be partly resolved with our MIDI baselines (V^ ~ 
0.4 for = 10 mas, B - 120 m, /I = lOyum), we also need to 
model the intensity distribution across the atmospheric layers. 
For these purposes, we use the dust-free dynamic model atmo- 
spheres, based on sel f -excited pulsation models (P and M se- 
ries), by ' Ireland et alJ (l2004bllaL and references therein) as the 
currently best available option. 

The P and M series were constructed to reproduce the M- 
type Mira prototypes o Cet and R Leo. The hypothetical nonpul- 
sating "parent" stars have masses M/Mq = 1.0 (P series), and 
1.2 (M series), Rosseland radius Rp/Re = 241 (P) and 260 (M), 
and effective temperatures T^g = 2860 K (P) and 2750 K (M). 
The models have solar abundances, luminosity L/Lq = 3470, 
and a pulsation period of 332 days, close to the periods of o Cet 
(332 days) and R Leo (310 days). They are available for 20 
phases in three cycles for the M series, and 25 phases in four 
cycles for the P series. These models have been shown to be 
consisten t with broad-band near -infrared interferom etric data of 
oCgt bv IWoodruff et alJ (l2004 and of R Leo by Fede le et all 
(l2005j) respectively. They have also been shown to be consis- 
tent with spectrally resolved, near-infrared interferometric data 
of the longer period (~ 430 days) Mira variable S Ori exhibit- 
ing a significant variation in the angular size as a function of 
wavelength, which was interpreted as the effect of molecular 
layers lying above th e continuum-forming la y ers as pr e dicted 
by these model series dWitfltowski et alJl2008l) . lTei et all (l2003h 
indicate that the spectral shape of Mira variables can be repro- 
duced reasonably well even for stars with different pulsation pe- 
riods. In comparison with the stellar parameters of the P and 
M model series, RR Aql pulsates with a longer period of 394 
days. T he main sequence pre cursor mass of RR Aql is 1.2 + 0.2 
M/Ma ('Wvatt&C ahnlll983l) . and its spectral type is M6e-M9 
(Samus et al. 2004), versus M5-M9 (P series) and M6-M9.5 (M 
series). For the purpose of the modeling effort conducted here, 
we can assume that the general properties of the P and M model 
series are also valid for a longer period variable such as RR Aql. 
However, an exact match of model parameters and observed pa- 
rameters as a function of phase may not be expected. 

The flux and visibility values of the combined atmosphere 
and dust shell ("global") model are computed for each spectral 
channel as 



/total — /star + /dust 



ratt 

V"* star 
total - 



ratt f 

-'star ^, Jdust,, 

"7 'star + "7 » dust 

/total /total 



(1) 

(2) 



where /!J"|. represents the attenuated flux from the dust-free 
model atmospherqfl (extended by a blackbody approximation for 
wavelengths beyond 23 //m), /dust represents the flux from the 
dust shell alone, /otai is the addition of these two terms, and Vstar 



' Note that /jJJJ,. may slightly vary across the stellar disk, which is not 
taken into account here. 
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Table 2. Best-fitting model parameters for each epoch 



Epoch Ovi 



(sil.) 



ii/-/^Phol 

(sil.) 



P 

(sil.) 



©Phot 

[mas] 



A 


0.58 


3.0 + 0.9 


4.5 ± 1.1 


2.5 + 0.4 


6.6 + 2.4 


B 


0.81 


3.0 + 0.7 


4.5 ± 1.6 


2.5 + 0.5 


7.1 +2.9 


C 


0.86 


2.5 ± 0.8 


5.0 ±1.3 


2.5 + 0.5 


8.0 + 2.8 


D 


1.53 


3.7 ± 0.6 


3.2 ± 1.1 


2.5 + 0.5 


6.7 + 2.6 


E 


1.77 


2.0 + 0.8 


4.5 + 1.4 


3.0 + 0.4 


8.0 + 2.9 


F 


2.46 


2.7 ± 0.8 


4.6+ 1.1 


3.0 + 0.5 


7.7 + 2.4 


G 


2.55 


3.0 ± 0.6 


3.4 ± 1.7 


2.1 +0.6 


7.2 + 2.5 


H 


2.62 


2.0 ± 0.8 


4.5 ± 1.1 


2.4 + 0.5 


8.5 + 2.7 


I 


2.74 


2.0 ± 0.9 


4.0 ± 1.5 


2.6 + 0.5 


8.5 + 2.8 


J 


2.84 


3.2 ± 0.7 


3.0 ± 1.5 


2.0 + 0.6 


7.2 + 2.6 


K 


3.43 


5.0 ± 0.8 


3.0 ± 1.3 


3.0 + 0.5 


7.3 + 2.2 


L 


3.54 


2.3 ± 0.9 


5.0+1.5 


2.5 + 0.4 


8.2 + 2.4 


M 


3.57 


2.5 + 0.7 


4.5+ 1.1 


2.5 + 0.4 


7.5 + 2.4 



Table 3. 


Average model parameters 






Model 


Tv Tv /^inZ-Z^Phot 

(AI2O3) (silicate) (silicate) 


P 

(silicate) 


©Phot 

[mas] 


M21n 


0.0 2.8 + 0.8 4.1+0.7 


2.6 + 0.3 


7.6 + 0.6 



and Vdust are the synthetic visibilities computed for the dust-free 
model atmosphere and for the dust shell, respectively. 

5. 1 . MIDI model parameters 

A clear silicate feature is identified in the spectra of RR Aql (see 
Fig.[3] and Sect. [Hi. RR Aql is one o f 31 ox ygen-rich stars stud- 
ied by iLorenz-Martins & Pompeial (l2000h . The authors mod- 
eled the stars using AI2O3 and silicate grains and suggest that 
the dust chemistry of RR Aql contains only silicate grains. In 
this stud y we examined both of these two dust s pecies, AI2O3 



grains (.Begeman n et al 
grains JOssenkop f et al. 



19971: lKoikeetal.l[T995l) and silicates 



1992). The amount of dust is described 



by the optical depths at Aq = 0.55 yum. The grain size was set 
to 0.1 ;um for all grains. We set the photospheric radius to the 
well-defined continuum photospheric radius at /I = 1.04/im 
(^Phot - ^1.04). The density gradient was defined by a single 
power law p(r) cc r^P with index p. The shell thickness was set 
to/?out/-R,., = 1000. 

The global model includes seven parameters. The radiative 
transfer model describing the circumstellar dust shell includes 
six parameters: the optical depths Ty (AI2O3) and Tv(silicate), 
the inner boundary radii RmlRp'aot (AI2O3) and RmlR^hoi (sili- 
cate), and the density gradients pA (AI2O3) and p^ (silicate). 
The model of the P/M atmosphere series is another parame- 
ter. Here, we used ten models covering one complete cycle of 
the M series: M16n (model visual phase Omodei = 0.60), M18 
(0.75), M18n (0.84), M19n (0.90), M20 (0.05), M21n (0.10), 
1VI22 (0.25), M23n (0.30), M24n (0.40), and M25n (0.50). 

We computed a grid of dust-shell models for each of these M 
models, including all combinations of optical depths Ty (AI2O3) 
= 0.0, 0.1, 0.2, 0.5, 0.8; Ty (silicate) = 2.0, 2.5, 3.0, 3.5, 4.0, 
4.5, 5.0; /?in//?Phot (AI2O3) = 2.0, 2.5, 3.0; RJRp^ot (silicate) = 
2.5, 3.5, 4.5, 5.5, 6.5; /?a (AI2O3) = 2.0, 2.5, 3.0, 3.5; and /^b 
(silicate) = 2.0, 2.5, 3.0, 3.5. 

In a first selection of suitable models, we compared the MIDI 
data of each epoch to our entire grid of models. The angular di- 
ameter ©Phot was the only free parameter. We increased the un- 
certainty of the part of the photometric spectra around 9.5 yum, 
which is strongly affected by telluric absorption. Otherwise, the 



weight of each data point was given by the corresponding uncer- 
tainty. For each epoch we kept the best ~30 models with lowest 
X^ values. In a next step, a grid with finer steps around these 
parameters was computed. This procedure was repeated several 
times. The results of the automatic selection were visually in- 
spected. We completed the selection with ten models for each 
epoch that agrees with the data best. 



5.2. MIDI model results 

For each of the 13 epochs we found the best-fitting model pa- 
rameters that include the best-fitting set of dust parameters of 
the radiative transfer model, together with the best-fitting model 
of the P/M series. We used the procedure outlined in Sect. 15. II 

Following iLorenz-Martins & Pompeial (l2000h . we investi- 
gated dust shells including AI2O3 and/or silicate shells with dif- 
ferent inner radii and density gradients. We obtained best-fit re- 
sults with a silicate shell alone, and the addition of an AI2O3 
shell did not result in any improvement in the model fits . This 
result is consistent with lLorenz-Martins & Pompeial (l2000l) . who 
classify RR Aql as a source that can be described with a silicate 
dust shell alone. 

The best-fitting parameters for each epoch are listed in 
Table |2| The table lists the epoch, the phase at the epoch, the 
optical depth Ty, the inner boundary radius /?in//?phot, the den- 
sity distribution p, and the continuum photospheric angular di- 
ameter ©Phot. Here, the dust shell parameters are those of the 
silicate dust shell. The errors of the dust shell parameters are 
derived in the same way as the standard deviation based on the 
find best-fitting models (Table|2]l. Table[3]lists the average model 
parameters of our different epochs. The average phase of our ob- 
servations is <i>v = 0.64 + 0.15. 

The best-fitting model atmosphere was M21n (r,nodei = 
2550 K, 'I>Modei=0.1) for all epochs covering minimum to pre- 
maximum pulsation phases (0.45-0.85). The difference between 
the average phase of our observation and the phase of the best- 
fitting model atmosphere can most likely be explained by the 
different stellar parameters of RR Aql compared to those of the 
M series. It means that Tf^g of RR Aql at the time of observation 
is similar to T^ff of model M21n, which has a different phase, but 
its basic parameters are also not exactly those of RR Aql. 

Figure |5] includes the model flux and model visibility com- 
pared to the observed values for the example of epoch A. For 
most epochs, the agreement between the models and the ob- 
served data is very good, in particular for the visibility spectra. 
The visibility spectra contain information about the radial struc- 
ture of the atmospheric molecular layers and the surrounding 
dust shells. In the wavelength range from 8/vm to ~9/im, the 
dust is fully resolved with baselines longer than ~ 20 m, and the 
visibility contribution of the dust reaches into the second lobe 
of the visibility function. RR Aql is characterized by a partially 
resolved stellar disk that includes atmospheric layers, with a typ- 
ical drop in the visibility function ~ IQjum, where the flux contri- 
bution of the silicate emission is highest and the flux contribution 
of the star relative to the total flux decreases. Beyond ~ 10/im, 
the spatially resolved radiation from the optically thin dust shell 
starts to be a notable part of the observed total flux. From ~ 
lOyum to 13/im, the dust contribution becomes nearly constant 
while the stellar contribution increases, and this results in a re- 
bound of the visibility function. Photometry spectra also fit well. 
For some epochs, the models predict higher fluxes than observed 
near the silicate emission feature at ~9.8jum. Generally, our at- 
tempt at a radiative transfer model of the circumstellar dust shell 
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that uses dynamical model atmosphere series to describe the cen- 
tral stellar source can reproduce the shape of both the visibility 
and the photometry spectra very well. Our model also follows 
the shape of the SED in the range of 1-40 /im (Fig.O. 

The obtained model parameters for the different epochs do 
not indicate any significant dependence on phase or cycle. This is 
consistent with the result from Sect. I4.1l that a direct comparison 
of visibility values of different phases or cycles did not show 
any variability within our uncertainties and that the photometric 
values indicated only small variations of up to ~ 2 cr. However, 
we cannot confirm or deny an eventual phase dependence of the 
dust formation process that affects the 8-1 3 //m visibility and 
photometry values by less than our uncertainties of about 5-20% 
and 10-50%, respectively. 

A silicate dust shell alone, i.e. without the addition of an 
AI2O3 dust shell, provides the best agreement with our data. The 
average optical depth of the silicate dust is Tv/(silicate)=2.8+0.8 
at /{ = 0.55 yum (corresponding to 0.03 at /I = 8//m, 0.06 
at /I = 12//m, and a maximum within 8-12//m of 0.22 at 
A - 9.8 /vm). The inner radius of the dust shell expressed in 
TJphot at /i = 1.04yum is Rm-AA+Q.! Rp^oi. The intensity profile of 
the dust-free dynamic model atm osphere extends to abou t 1 .5- 
2 photosphere radii at lOjum (cf. 'Wittko wski et al.|[2007h . The 
power-law index of the density distribution is p=2.6+0.3. Wind 
models predict a power-law index of the density distribution of 
2 at radii outward the dust formation zone. However, for radii 
close to the dust formation z one (r<~10 R*), s. larg er index is 
expected (cf, e.g.. Fig. 6 in IWittkowski et alJf2007h . which is 
consistent with our result. 

The average photospheric angular diameter results in ^hot - 
7.6 ± 0.6 mas. This value is lower than the K-band {A - 2.2//m, 
^A = 0. 4 ^m) UP diarneter of ©ud = 10.73 + 0.66 mas de- 
rived by Ivan Belle et alj (|2002|) at a minimum phase of 0.48. 
This can most likely be explained by the different radius def- 
initions. Our diameter is a photospheric angular diameter that 
has already been corrected for effects of molecular layers lying 
above the continuum photosphere using the prediction by the 
best-fitting model atmosphere, while the diameter by van Belle 
et al. is a broad-band uniform disk diameter that still contains 
contamination by molecular layers. It is known that broad-band 
UD diameters may widely overestimate the photospheric radius 
owing to contaminati on by molecular lay e rs lying above the 
photosphere (cf. e.g Feland et al.l l2004bt iFedele etall l2005t 
IWittk owski et al.'2007^. With the parallax of tt = 1.58 +0.40 mas 
dVle mmings & van Langevelde 2007), our value for 0phot corre- 
sponds to a photospheric radius of /?phot = 52O^j4g7?0. Together 
with the bolometric ma gnitude of RR Aql of nzboi = 3.71 
and Amboi = 117 from IWhitelock et alj ( l2000l) . our value for 
^hot at Ov ~ 0.64 corresponds to an effective temperature of 
Jeff ~ 2420 ± 200 K. This value is consistent with the effective 
temperature of the best- fitting model atmosphere M21n, which 
isTeff = 2550K. 



5.3. Model simulations 

We performed a number of model simulations in order to inves- 
tigate the visibility and photometry variations that are theoret- 
ically expected in the 8-1 3 //m wavelength range for the typi- 
cal parameters of RR Aql as determined above in Sect. 15.21 We 
based our simulations on a global model with typical parameters 
similar to those derived for RR Aql (cf. TableO, and investigated 
the effects of expected variations in the atmosphere model and 
the dust shell parameters during a pulsation cycle on the observ- 
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Fig. 10. Simulation 1 . Synthetic flux (top) and visibility ampli- 
tude (bottom) in the wavelength range 8-13jum. The solid lines 
represent the global models, the dashed-dotted lines denote the 
unattenuated stellar contribution i.e. Vstai in Eq. |2] the dotted 
lines denote the attenuated stellar contribution, and the dashed 
lines denote the dust shell contributions. This simulation com- 
pares two models consisting of the same dust shell parameters 
but different model atmospheres, the post-maximum model at- 
mosphere M2 In (i^niodei - 0.10), and the minimum model atmo- 
sphere M23n (Omodei - 0.30). The photospheric angular diame- 
ter is assumed to be larger at post-maximum phase than at mini- 
mum phase (8.5 mas compared to 7 mas). The projected baseline 
length is 60 m. For the exact model parameters, see Table|4] They 
describe variations around our best-fitting model for RR Aql. 



able photometry and visibility spectra. We mainly used a mean 
projected baseline length of 60 m. As in Sect. |5] the global model 
consists of a radiative transfer models of the AI2O3 and silicate 
dust shells where the central source is described by a dust-free 
dynamic model atmosphere. Table |4] lists the performed simu- 
lations including the phase of the model, the optical depth ry 
(AI2O3) and Ty (silicate), the inner boundary radii of the dust 
shell ^in/^phot (AI2O3) and /?in//?phot (silicate), the power-law 
indices of the density distributions p (AI2O3) and p (silicate), 
and the continuum photospheric angular diameter ©phot- We as- 
sume a phase depende nce of the angular photospheric diame- 
ter of about 20% (Tho mpson eta l. 2002a; Ireland et al. 2004b|). 
Figures [T0l-[T4lshow the results from the simulations. Each sim- 
ulation compared two models, where one or more of the model 
parameters have been varied. 

Simulation 1 (Fig. [TOl i In simulation 1, we investigate the effect 
of model atmospheres at different pulsation phases and effective 
temperatures. We compared model M21n at post-maximum pul- 
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Table 4. Model simulations, comparing two models that differ in one or more parameters (marked by bold face). 



Simulation 


Model 


Ovis 


Tv 

AKjOs) 
tva 


Tv 

(silicate) 
tvb 


(AI2O3) 

ra 


RmlRvhol 

(silicate) 
rb 


P 

(AI2O3) 

pa 


P 

(silicate) 

pb 


©Phot 

[mas] 
diam. 


Projected 
baseline [m] 

Bn 


1 


M21n 
M23n 


0.1 
0.3 


0.0 
0.0 


3.0 
3.0 


2.2 
2.2 


4.5 
4.5 


2.5 
2.5 


2.5 
2.5 


8.5 

7.0 


60 
60 


2 


M21n 
M23n 


0.1 
0.3 


0.0 
0.0 


3.0 
4.0 


2.2 
2.2 


4.5 
3.8 


2.5 
2.5 


2.5 
2.5 


8.5 
7.0 


60 
60 


3a 


M21n 
M21n 


0.1 
0.1 


0.0 
0.0 


3.0 
5.0 


2.2 
2.2 


4.5 
3.0 


2.5 
2.5 


2.5 
2.5 


8.5 
8.5 


60 
60 


3b 


M21n 
M21n 


0.1 
0.1 


0.0 
0.0 


3.0 
5.0 


2.2 
2.2 


4.5 
3.0 


2.5 
2.5 


2.5 
2.5 


8.5 
8.5 


40 
40 


3c 


M21n 
M21n 


0.1 
0.1 


0.0 
0.0 


3.0 
5.0 


2.2 
2.2 


4.5 
3.0 


2.5 
2.5 


2.5 
2.5 


8.5 
8.5 


25 
25 


3d 


M21n 
M21n 


0.1 
0.1 


0.0 
0.0 


3.0 
5.0 


2.2 
2.2 


4.5 
3.0 


2.5 
2.5 


2.5 
2.5 


8.5 
8.5 


15 
15 


4 


M21n 
M21n 


0.1 
0.1 


0.0 
0.6 


3.0 
3.0 


2.2 
2.2 


4.5 
4.5 


2.5 
2.5 


2.5 
2.5 


8.5 
8.5 


60 
60 



Notes. Here, Ovjs is the phase, Ty the optical depths, i?m//?phot the inner boundary radii, p the density gradients, Qphot the angular diameter. 
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Fig. 11. Simulation 2. Like simulation 1 in Fig.[TOl but for the 
parameters of simulation 2 comparing a post-maximum atmo- 
sphere model and a minimum model as in simulation 1, but 
where the dust is also assumed to be closer to the star with larger 
optical depth at minimum phase and farther from the star with 
lower optical depth at post-maximum phase. 



sation phase (Omodei - 0.10, Teff = 2550 K) and model M23n at 
pre-minimum pulsation phase ((tmodei = 0.30, T^ff = 2230 K), 
applying the same dust shell parameters for both models. Here, 
the phase difference between M21n and M23n of 0.2 corre- 



Fig. 12. Simulation 3. As for simulation 2 in Fig. \TT\ but where 
the atmosphere model is not varied. 



sponds to the phase difference of our observations (-0.64+0.2) 
The dust parameters were based on the average parameters (see 
Tab. O of RR Aql derived in Sect. 15.21 The projected baseline 



length was set to Bp = 60 m, which is a typical mean value 
of our MIDI observations. The overall shapes of the visibil- 
ity spectra Vtotai (M21n) and Vtotai (M23n) are similar. The dif- 
ferences between the models are wavelength-dependent in the 
range of 5-20%. The shape of the photometry spectra are sim- 
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Fig. 13. Simulation 3. As Fig. [12] but showing the visibihty spec- 
tra for different projected basehne lengths of (from to bottom) 
40 m, 25 m, and 15 m. 



Fig. 14. Simulation 4. Comparison of a model with a silicate 
shell only and a model that also includes an AI2O3 dust shell 
with lower optical depth. For the model parameters, see Tabled 
For the description of the figure, see the caption of Fig. [TO] 



cated more dust formation near minimum pulsation phase, with 
inner boundary radii of the dust shell located closer to the star. 
We increased the optical depth for model M23n. We also set a 
more extended dust shell after visual maximum for model M21n 
in comparison to the model M23n. This setting results in spec- 
tra with differences in the range of 2-19% for the visibility val- 
ues. The photometry differs again mostly around lO/zm with a 
maximum difference of 20%. Compared to simulation 1, the re- 
sults are very similar, but with less differences in the photometry 
spectra. As for simulation 1, the results are consistent with our 
not detecting visibility variations and marginally detecting pho- 
tometry variations within our measurement uncertainties. 



ilar for both tested models /total (M21n) and /total (M23n). The 
largest difference of ~ 25% is seen around 10 //m, with less flux 
at minimum pulsation phase. Compared to the uncertainties of 
our MIDI measurements of ~ 5-20% in the visibility spectra and 
of ~ 20% on average in the photometric spectra, simulation 1 is 
consistent with our non-detection of visibility variations and a 
marginal detection of photometry variations. 

Simulation 2 (Fig. [TTI) Compared to simulation 1, we adjusted 
the parameters according to the assumption that the outer lay- 
ers of the atmosphere are cooler near minimum visual pulsa- 
tion phase, and therefore more dust grains can form, and a 
higher molecular opacity can be expected. This theoretical as- 
sumptio n is consistent with mid -infrared interferometric obser- 
vations (IWittkowski et al.ll2007h . where the observed data indi- 



Simulation 3 (Figs.[T2]and[T3ll Here, we investigate the effect of 
different dust shell parameters (as in simulation 2), but keeping 
the M model constant. We used model M21n (phase 0.1), which 
was the best-fitting model to our RR Aql epochs and also used a 
constant photospheric angular diameter. We compared two (sili- 
cate) dust models as in simulation 2, where in one model the op- 
tical depth is lower and the inner radius of the dust shell larger 
compared to the other model. These sets of model parameters 
lead to very similar results, where the visibility spectra almost 
overlap, with a maximum difference of 4% for a projected base- 
line length of 60 m (Fig. fT2l i. The photometry spectra are also 
very similar with a maximum difference of 7%. Figure[T3]shows 
the visibility results for different projected baseline lengths of 
40 m, 25 m, and 15 m. With projected baseline lengths of less 
than 60 m, the difference between the visibility values become 
more pronounced (10% difference beyond lOyum with projected 
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baseline lengths of 40 m and 25 m). This can be understood be- 
cause the dust shell is over-resolved with a projected baseline 
length of 60 m, so that the visibility cannot be sensitive to vari- 
ations in the geometry of the dust shell. With lower projected 
baseline length, the visibility values are located in the first lobe 
of the visibility function, hence more sensitive to variations of 
the dust shell geometry. However, with the even lower projected 
baseline lengths of 15 m, the difference between individual mod- 
els again starts to be less pronounced (<7%). This indicates an 
optimum projected baseline length of ~ 25^0 m to characterize 
the dust shell geometry for sizes as found for RR Aql. 

Simulation 4 (Fig.[T4|i Furthermore, we investigate the effect of 
adding small amounts of AI2O3 grains to a thus far pure silicate 
dust shell. We examined Ty(Al203) of 0.3 and 0.6. Figure [14] 
shows the simulation with Tv(Al203) of 0.6, showing that the 
addition of small amounts of AI2O3 grains leads to almost iden- 
tical photometry and visibility spectra (differences <2%). 



6. Discussion 

Our MIDI observations of RR Aql did not show significant vari- 
ations in the 8-1 3 /im visibility within the examined pulsation 
phases between 0.45 and 0.85 within a total of four pulsation cy- 
cles, and showed only marginal variations in the 8-1 3 /im flux. 

We performed model simulations with expected variations 
of the pulsation phase of the innermost dust-free atmosphere 
model and of the parameters of the surrounding dust shells. 
These model simulations show that visibility variations are in- 
deed not expected for the parameters and observational settings 
of RR Aql at wavelengths of 8-1 3 //m within the uncertainties 
of our observations. Variations in the flux spectra may in some 
cases just be detectable. Thus, our observational result of a con- 
stant visibility and only slightly varying flux at wavelengths of 
8-13jum are consistent with, and not contradicting, theoretical 
expectations of a pulsating atmosphere. 

Our model simulations indicate that detections of pulsation 
effects at mid-infrared wavelengths would in particular benefit 
from smaller uncertainties in the photometric spectrum than in 
our current data. Also, a wide range of projected baselines at 
each phase, for RR Aql in particular including baseline lengths 
around 20-30 m, would help us to distinguish models of different 
pulsation phases. 

For our analysis, we have considered silicate (Ossenkopf et 
al. 1992) and AI2O3 (" Beeemann etap '1997^. 'Koike et al."l995') 
dust species, following the work byJLorenz-Martins & Pompeia 
(I2OO O). and used fixed grain sizes of 0. 1 /im. It is known that 
other dust species with a more complex gran size distribu- 
tions may occur in the circumstellar environmen t of Mira vari- 
ables (e.g.. lHofned2008l:lMolster & Watersl2003l and references 
therein). We showed that a model including only a silicate dust 
shell can reproduce the observed RR Aql visibility and flux spec- 
trum at 8-1 3 /im. The addition of an AI2O3 dust shell with com- 
parable low optical depth did not significantly improve the fit to 
our data. However, our model simulations have shown that our 
8-1 3 //m visibility and flux values are not sensitive to the addi- 
tion of an AI2O3 dust shell with low optical depth within our un- 
certainties. As a result, we cannot exclude the presence of an in- 
ner A I2O3 dust shell in addition to the silicate dust shell. IWoitkd 
(1200 6) uses dynamical models for dust-driven winds of oxygen- 
rich AGB stars, including frequency-dependent radiative trans- 
fer, and finds that dust temperatures strongly depend on material. 
Two dust layers are formed in his dynamical models, almost pure 



glassy AI2O3 close to the star (r > 1 .5 ^star) and the more opaque 
Fe-poor Mg-Fe-silicates farther out at 4-5 Rstai- 

The dust-free dynamic model atmospheres predict a signif- 
icant dependence of the characteristics of the molecular lay- 
ers on the stellar puls ation phase at near-infrared wavelengths 
([Ireland et al.ll2004blah . There are a few near-infrared interfero- 
metric observations, which detected a clear variation in the con- 



tinuum a ngular diameters with pulsation phase (Thompson et al. 



2002b; Perrin et all Fl999l lYoung et all l2000bt IWoodruff^et al. 



2004) . Spectrally resolved near-infrared interferometric mea- 
su rements at different phas es, such as the AMBER observations 
bv Wittk owski et all (|2008") but at more than one phase, promise 
to lead to stronger constraints on dynamic model atmospheres. 

Our model simulations based on the combination of dy- 
namic model atmospheres with a radiative transfer model of 
the dust shell predict only small variations with phase at mid- 
infrared wavelengths. The difference near a wavelength of 10/im 
amounts to ~5% for the visibility values and to ~25% for the 
photometry values. These results lead to the suggestion that 
the stellar photosphere and overlying molecular layers pulsate, 
which is demonstrated by the diameter variations in the near- 
infrared, but that these pulsations cannot be detected for RR Aql 
by our 8-1 3 /im interferometry within our uncertainties. In addi- 
tion, our visibility uncertainties do not allow us to exclude vari- 
ations in the geometry and optical depth of the dust shell as a 
funct ion of pulsation ph ase as observed previously for other tar- 
gets (iLopez et al.llI997h . Possible explanations for not detecting 
a phase-dependence of the dust shell parameters in our study 
compared to the ISI observations include: (i) longer baseline 
lengths in our study that often fully resolve the very extended 
silicate dust shell; (ii) our limited phase coverage between mini- 
mum and pre-maximum (0.45-0.85) phases. Possible variations 
over the whole pulsation cycle cannot be excluded. 

The observed variability of the 8-13 //m flux at a significance 
level of l-2cr may indicate variations in the stellar radiation re- 
emitted by the dust and/or changes in either the dust geometry 
or optical depth. 

7. Summary and conclusions 

We have investigated the circumstellar dust shell and character- 
istics of the atmosphere of the oxygen-rich Mira variable RR Aql 
using mid-infrared interferometric observations. We observed 
RR Aql with the VLTI/MIDI instrument at different pulsation 
phases in order to monitor the photometry and visibility spectra. 
A total of 57 observations were combined into 13 epochs cover- 
ing four pulsation cycles between April 2004 and July 2007, and 
covering pulsation phases between minimum and pre-maximum 
phases (0.45-0.85). 

We modeled the observed data with an ad-hoc radiative 
transfer model of the dust shell usi ng the radiative transfer code 
mcsimjnpi bv lOhnaka et al.l(l2006h . In this way, we used a series 
of dust-free dynamic model atmospheres based o n self-excited 
pulsation models (M series. [Ireland et aDl2004bllah to describe 
the intensity profile of the central source. This study represents 
the first comparison between interferometric observations and 
theoretical models over an extended range of pulsation phases 
covering several cycles. 

Our main observational results are as follows 

- The interferometric data do not show any evidence of intra- 
cycle visibility variations. 

- The data do not show any evidence of cycle-to-cycle visibil- 
ity variations. 



13 



Karovicova et al.: Mid-infrared interferometric monitoring of the Mira variable RR Aql 



- The 8-13|/m flux suggests intracycle and cycle-to-cycle 
photometry variations at a significance level of l-2cr. 
Follow-up observations with higher accuracy using a dedi- 
cated photometric instrument, such as VISIR at the VLT, are 
needed to confirm this result. 

These observational results can be explained by dynamic 
model atmospheres and variations in the dust shell parameters. 
Simulations using different phases of the dynamic model atmo- 
sphere and different sets of dust shell parameters predict vis- 
ibility variations that are lower than or close to the observed 
visibility uncertainties (5-20 %). Model-predicted variations of 
the photometry spectra are largest around wavelengths of lOjjm 
with difference of up to ~25% corresponding to up to l-2cr. 

The best-fitting model for our average pulsation phase of 
Oy - 0.64 + 0.15 includes a silicate dust shell with an op- 
tical depth of Tv(silicate) - 2.8 + 0.8, an inner radius of 
Rin = 4.1 + 0.7 Rphoi, and a power-law index of the density dis- 
tribution of p = 2.6 + 0.3. The corresponding best-fitting at- 
mosphere model of the series used to describe the central in- 
tensity profile is M21n (Tj^odei - 2550 K, <E)modei - 0.1) with 
a photospheric angular diameter of flphot = 7.6 + 0.6 mas. The 
photospheric angular diameter corresponds to a photospheric 
radius of Rphot = 52Q^j^^Rq and an effective temperature of 
Teff ~ 2420+200 K. The latter value is consistent with the effec- 
tive temperature of the used model M21n. The combined model 
can reproduce the shape and features of the observed photome- 
try and visibility spectra of RR Aql well, as well as the SED at 

We conclude that our RR Aql data can be described by a sil- 
icate dust shell surrounding a pulsa ting atmosphere, consistent 
with observations of Mira variables jLorenz-Martins & Pompeial 
I2OOQ1) . The effects of the pulsation on the mid-infrared flux and 
visibility values are expected to be less than about 25% and 20%, 
respectively, and are too low to be detected within our measure- 
ment uncertainties. Although the addition of an AI2O3 dust shell 
did not improve the model fit, our simulations also indicate that 
we cannot exclude the presence of an inner AI2O3 dust shell with 
relatively low optical depth, which may be an important contrib- 
utor to the dust condensation sequence. 

Our modeling attempt at a radiative transfer model of the 
dust shell surrounding a dynamic dust-free model atmosphere 
provides constraints on the geometric extensions of atmospheric 
molecular and dust shells. It should be noted that it cannot 
explain the mechanism by which the observed mass loss is 
produced and the wind is driven. Newer models of M type, 
i.e. of oxygen-rich Miras, that include atmospheric dust (e.g. 
[Ireland et al.ll2008l Ireland, in preparation) typically predict ef- 
fective acceleration very close to zero in high layers but, so far, 
no clear outward acceleration, i.e. no wind. Subtle details are be- 
ing discussed that may over come this sho rtcoming of dynamic 
models of M type Miras (e.g. lHofnerl201 ih . Future observations 
aiming at characterizing and constraining other new models of 
the mass-loss process and the wind driving mechanism at mid- 
infrared wavelengths would benefit from obtaining more precise 
photometry values using a dedicated instrument like VISIR at 
the VLT, the addition of shorter baselines to characterize the ex- 
tension of the silicate dust shell, and a more complete coverage 
of the pulsation cycle. The addition of concurrent, spectrally re- 
solved, near-infrared interferometry would be needed to more 
strongly constrain atmospheric molecular layers located close to 
the photosphere. 
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Fig. 5. VLTI/MIDI interferometry at 8-13//m of RR Aql for the example of epoch A (stellar phase 0.58) and epoch I (stellar phase 
2.75) . The panels show (a) the flux, (b) the visibility amplitude, (c) the corresponding UD diameter, (d) the corresponding Gaussian 
FWHM diameter as a function of wavelength. The gray shading indicates the wavelength region around 9.5 //m that is affected by 
atmospheric absorption. Panel (e) shows the visibility amplitude as a function of spatial frequency for three averaged bandpasses of 
8-9 ;um, 10-1 1/im, and 12-13/im. The crosses with error bars denote the measured values. The solid lines indicate our best-fitting 
model, as described in Sect. |5] It consists of a combination of a dust-free dynamic model atmosphere representing the central star 
and a radiative transfer model representing the surrounding dust shell. The contributions of the stellar and dust components alone 
are indicated by the dotted and the dashed fines, respectively. 



